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Abstract 

The average thermodynamic power of a time-dependent external potential in the white-noise 
Langevin model is derived using a Green's function solution. The power appears as a driving term 
in the differential equation for the average energy and determines whether the solution is stationary 
or non-stationary. Different dynamics are illustrated with explicit models: a linear potential with 
a static magnetic field, a linear potential perturbed with an oscillating component and a magnetic 
field switch modeled using a tanh protocol. 
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I. INTRODUCTION 



Time-dependent potentials in the Langevin equation have been thoroughly used by Jayan- 
navar et al. in investigating properties of non-equilibrium systems, a few examples of which 
are Refs. In these papers the distribution of thermodynamic work done on the 

Brownian particle by the external potential is analytically or numerically calculated and is 
then analyzed using non-equilibrium criteria such as the Jarzynski and Crook relations. The 
time-dependence of the potential facilitates the reversal of the protocol, which is central in 
comparing forward and backward work statistics. Moreover, the time-dependence sets the 
study of non-equilibrium statistical mechanics in the necessary context of non-equilibrium. 
As we will see, time-dependence is essential in producing non-stationary solutions. 

Aside from the context of non-equilibrium theorems, understanding time-dependence is 
also interesting from a phenomenological point of view. In particular, knowledge of heavy 
ion and plasma thermodynamics in the presence of a time-dependent electromagnetic field 
could lead to methods of experimental control. An example is the controlled diffusion of 
overdamped Brownian particles using a time-dependent drift term {4]. 

In this work we investigate the role of a time-dependent potential ~ g{t)r'^ in the white- 
noise Langevin model by explicitly evaluating the mean and variance of the phase space 
variables. We derive the formal solution using the Green's function method and pay partic- 
ular attention to a differential equation describing the time evolution of the average energy. 
The dynamics is driven by the average thermodynamic power, which is found to be propor- 
tional to the product the time-dependent potential strength and the diffusion coefficient. 
Then it is the behavior of these two quantities that directly implies whether the system 
would be stationary or non-stationary in the long time limit. We discuss this and other 
physical implications by evaluating sample systems. In Sect. [IV]we briefly revisit a static 
system previously treated in Ref . [g] . We use an alternative approach where a direct connec- 
tion between the thermodynamic power and the mean behavior of the particle is established. 
In Sect JV] we consider a linear potential perturbed with an oscillating part and we find that 
the calculated thermodynamic quantities exhibit steady-state oscillations as well. In Sect. 
IVII we treat a time-dependent magnetic field and focus on a magnetic field switch using a 
tanh protocol. The system is evaluated by matching perturbation solutions from the ex- 
treme values of an expansion parameter. We find that a system in equilibrium is taken out 
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of equilibrium by the non-zero power during the time of switching and that the work done 
by the induced electric field is not the same for forward and backward switch protocols. 



II. LANGEVIN MODEL 

We consider a charged Brownian particle in the presence of a static magnetic field B = Bz 
and an external time-dependent potential ^g{t)r'^. The Langevin equation in the (x, y) plane 
can be written in terms of the matrix cr = ( i'l o ) as 

VV yy) yyj W yyj 

where u = qB/m is the cyclotron frequency for a particle of charge q and mass m, (3 is the 
damping factor and fj = rj^x + riyy is the stochastic force per unit mass. The Cartesian com- 
ponents of the stochastic force are Gaussian, have zero mean and obey the autocorrelation 

(r^,(tH.(t'))=5.,^5(t-t'). (2) 

If we rewrite Eq.([T]) using the eigenvectors of cr as a basis we can reduce our working 
equation to one-dimension: 

vi + {(3 - iuj]vi + g{t)ri = 7]i. (3) 

Recovering the real valued quantities can then be done using the definition qi = qx — iqy for 
q e {v,r,r]}. 



III. GREEN'S FUNCTION SOLUTION 



We solve Eq.Q by treating the stochastic force as an inhomogeneous and finding the 
Green's function (gI. The Green's function has the properties 



G{t,T) 



= and -^G{t,T] 
T=t dr 



-1 



T=t 



(4) 



and obey the equation of motion 



-^G{t, r) + {/3 - ^oJ}^G{t, r) + g{r)G{t, r) = 0. 



(5) 
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From Eqs. (jlj) and ([5]) the Green's function can be determined as 

Git,r)=mm^jm^, (6) 

/2(i)/iW-A(i)/iW 
where /i(t) and /2(t) are the hnearly independent solutions of Eq.Q. 

Then the solution can be written as 

ri = (ri)+ f G{t,T)r,{r)dT, (7) 







where the ensemble average (ri) is given by 



(ri) = bi(0) + ri(0) (/3 - lu)] G{t, 0) - n(0)|:G(t, r) . (8) 



r=0 



The velocity is just the derivative of Eq.([7]) 

vi = (vi) + [ H{t,T)7]{r)dT, (9) 
Jo 

where H{t,T) = iG{t,T). 

With the solutions written in this form we can use Eq.([2]) to arrive at the variance in 
position and velocity for i = x,y given respectively as 

{Ar^) = ^ fdr\Git,r)\^ (10) 
fn Jo 



2, 2kTp 



and 

{Avt) = ^-^ I dr\H{t,T){\ (11) 

where |/| denotes the modulus of /. And we can readily compute the diffusion coefficient 
from the definition -Dj(t) = ^{Arf) as 

A(t) = ^ /* Re [H{t, T)G*{t, r)] dr, (12) 
m Jq 

where /* denotes the complex conjugate of /. 

A differential equation for the time evolution of the average energy per unit mass (E). = 
\ [("^i^) + (^"^i^)] can be derived by multiplying Eq.([3]) by vl and then taking the thermal 
average, giving 



k 

(13) 
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In the derivation we used the property H(t,t) = 1, which can be deduced from Eq.Q. 
Note that there is no magnetic field contribution in Eg. (1131) since it cannot do work. The 
subscript k denotes the parts of the equation which carry a factor of the Boltzmann constant 
k ~ 10^^^. Hence Eq.f lT^ can be seen as two equations separated by this scale. 
The macroscopic equation is 

~{v^f + Piv.f = -^-^Re [{n){vl)] . (14) 

The term on the right hand side is the average power of the external force. Without this 
term the kinetic energy would vanish exponentially as exp[— 2/5t] due to damping. 
The microscopic equation, 

l^(Av^) + f3{Avl) - ^ = f Re[H{t,r)G%t,r)]dr (15) 

is similar except for the constant term — on the left hand side. The term on the right 
hand side is the average thermodynamic power. Without this term the variance in velocity 
would decay as exp[— 2/3t] to (Aff ) = ^ and the average energy per unit mass would become 
J2i = kT/m. Note that the power can be written as P{t) = ^g{t)Di(t). This means 
that if the system were to equilibrate then either g{t) or Di{t) should vanish at large times. 
If both the potential and the diffusion coefficient were non-vanishing then the variance in 
velocity would not be stationary in the long time limit. The behavior depends on the type 
of driving that P(t) provides. 

When the work done is performed adiabatically such that the power is negligible, we are 
led to a relation derived by Mazo under the assumption of an equilibrium process [7|: 

Clearly this does not hold in the non-adiabatic case since the work done allows for an increase 
in kinetic energy (wj(t))^ > fi(0)^, which would make the variance negative. 

In the following sections we look at some examples illustrating the role of the average 
thermodynamic power in the statistics of a Brownian particle. 
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IV. STATIC LINEAR POTENTIAL WITH A STATIC MAGNETIC FIELD 



When g{t) = g the Green's function can be solved as 



G{t, r) = ^ (17) 

/i 



where fi = \/ —^g + (/5 — icu)^. With the Green's function at hand computing the thermo- 
dynamic quantities can be done following the previous analysis. However we could approach 
the system from a different perspective. By considering a change of variable in Eq.Q 
r ^ t' — s, the equation of motion for the Green's function along a time slice t = t' becomes 

t' - s) + {/3 - tiu}-^G{t', t'-s)+ gG{t\ t' - s) = 0. (18) 

Since s runs from to t' it follows from the properties of G{t, r) that the initial conditions 

are G{t', t' - s) =0 and £G{t', t' - s) = 1. 

s=0 s=0 

Note from Eq.Q that Eq. ffTSl) is now identical in form to the equation of motion for 
the average position. Further more, since we are interested in the long-time behavior we 
can choose the initial conditions (ri(0), f i(0)) = (0, 1) which correspond to G(0, 0) = and 
if (0,0) = 1. This reduces the average position given in Eq.(l8]) to (ri) = G(s,0), which 
obeys the equation of motion 

-^Gis, 0) + {/3- ^^}^G(s, 0) + gis)Gis, 0) = 0. (19) 

Since the equations of motion and the initial conditions of G(t', t — s) and G{s, 0) are the 
same, this establishes that G{s, 0) = G{t',t' — s) for s G [0,t']. Hence we find that 

[ Re[G*{t',T)H{t',T)]dT = I Re[G*{t',t' -s)H{t',t' -s)]ds 
Jo Jo 



t' 



Re [G*{s, 0)H{s, 0)]ds 



'0 

= i|G(t',0)p. (20) 

The left hand side of the above equation is proportional to the diffusion coefficient and hence 
power. While the right hand side is the square of the mean position of the particle. Thus 
when the mean position goes to zero in the long time limit, the diffusion coefficient and power 



also go to zero. Consequently the variance in velocity goes to kT/m according to Eq. (fT5!) . 
The exact behavior can be readily studied from the explicit Green's function given in Eq. ffTTl) . 
Note that the uniqueness of the phase space equilibrium point ((ri), (f i)) = (0, 0) guarantees 
that the equipartition value is the only stationary value for the velocity variance. 

If we evaluate Eq. fllOl) using Eq. fll7p and take the long time limit, we find that the variance 
in position approaches 

kT 

lim (Arh = —. (21) 

This defines a radius ro = ^J^^ within which the Brownian particle is most likely to be 
found. This result is independent of the static magnetic field, which only affects the rate 
at which the system diffuses to tq. A stronger magnetic field results in a slower diffusion 
process since a spiraling path results in less displacement. 

When g < the potential is repulsive and the mean becomes unbounded in phase space. 
Hence the diffusion coefficient and variance in velocity also increase indefinitely. 

V. OSCILLATING LINEAR POTENTIAL 

We consider a system with no magnetic field u = and an oscillating linear potential 
g{t) = g + ecos[Qt]. The solutions to Eq.(l5]) are 

4.-/?^ -2e gr- ^^^^ 



/i(r) = e^'/^MathieuC 
and 

Mr) = e^^/^MathieuS [^^, ^, f ] . (23) 

The Green's function is then constructed using Eq.(l6]) and when e = is set, we recover the 
Green's function given in Eq.( |T71) with u = 0. 

If e/g « 1 we can consider the oscillatory part to be a perturbation and take Eq. fl2T|) 
to lowest order in e as 

lim (Arh = — (l-- cosffitl^ . (24) 

This describes how the radius tq oscillates as it reacts to the instantaneous value of g{t). 
Since the diffusion coefficient does not go to zero i.e. -Dj(t) ~ nsin[r2t], the thermodynamic 



power P{t) ~ g{t)Di{t) also remains oscillating. Then Eq. lfTSi) implies that the variance 
in velocity oscillates about kT/m in the long time limit. In Figured] we show the numer- 
ically integrated thermodynamic quantities. In this example the harmonic strength g{t) is 
always positive and the mean position becomes stationary at zero at large times, while the 
diffusion coefficient oscillates. In this case the connection between the mean dynamics and 
thermodynamics derived in Eq. (!20|) does not hold. 



VI. MAGNETIC FIELD SWITCH 

In the presence of a time-dependent magnetic field B = B{t)z, an electric field of the 
form E = — f is induced. We can neglect further electromagnetic inductions because 

of the ~ 1/c scaling. The Langevin equation we are interested in solving is 

+ {(3 - tuj{t)}vi - tX'^n = r]i, (25) 

where A is a book-keeping parameter we will set to 1. For some choices of uj{t) the Green's 
function for the above equation can be solved using the previous formal solution. However, 
in this example we are interested in using the tanh function as a magnetic field switch 
protocol. We will solve for the Green's function by perturbation through constructing a 
matched series from the known A = and A = 2 solutions. 
From Eq. fl25|) the velocity vi can be expressed implicitly as 

v,=v,{0)To{t,0)+ [ ro(t,r)r/i(r)rfr + 4 / To{t,T)dj{T)n{r)dT, (26) 
Jo ^ Jo 

where we have the velocity Green's function in the A = case 

ro(t, r) = exp (- ^ [/^ - M^)] ■ (27) 
To express vi explicitly as in Eq.Q we iterate Eq.( l26i) . arriving at the Green's function 



p=0 

where 



rp{t,T)= [ w{t,s)rp_i{s,T)ds (29) 

Jo 
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and 



W{t,T) = / To{t,s)ilj{s)ds. 



(30) 

Similarly we can construct a perturbative solution about A = 2. Note that Eg. (1251) can 
be written as 



(31) 



Thus f 1 can be expressed implicitly as 

Vi=Vi{0)Ao{t,0)+ [ Ao{t,T)r]i{T)dT + l^^ [ Ao{t,T),dj{T)n{T)dT, (32) 

Jo ^ Jo 

where we have the Green's function for the A = 2 case 



d /•* 

Ao{t,r) = - ro{t,s)ds. 



Again by iteration we have the Green's function 



p=0 



where 



Ap(t,r) = / X{t,s)Ap.i{s,T)ds 
Jo 



and 



X{t,T) = j Ao{t,s)uj{s)ds. 
By combining Eqs. fl28|) and fl34l) we connect the two solutions as 



p=0 



ro(t,r)+Ao(t,r) 



We will use the zeroth order term 

as the velocity Green's function and 

G'o(t,r) = / HQ{s,T)ds 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



as the position Green's function. To check how well the approximation performs we compare 
the calculations for the field u!{t) = ,^tanh[t] with those of the exact solution given by the 
functions 



/i(r) = 2F1 

and 

/2(r) = 2F1 



1 + e 



2r 



(40) 



-1 



1 + e 



2t 



2 



2 ' 2 ' 2 ' 1 + e2* 

(41) 

where 2-F1 is a hypergeometric function. Figure [2] shows good agreement between the results. 
Note that the diffusion coefficients exhibitjiransient oscillations and their final values agree 
with the static magnetic field expression 



T^static _ ^^-^ fAcy\ 

~ m (3^ + 00^- ^^^^ 
The average thermodynamic power of the induced electric field is given by 

P{t) = Re[iu{t)D{t)], (43) 

where uj{t) = ^sech^[t]. A larger switch magnitude ^ induces a stronger electric field, result- 
ing to a larger transient spike in the velocity variance. 

Note however that the previous computations do not represent a well prepared system 
since the switching starts at t = when the transient exponential decay is still in effect. 
What we would like to see is switching in a system already in equilibrium as modeled by the 
protocol uj{t) = tuo + itaiih.[a{t — ts)] for ts larger than the exponential decay time In 
this case the diffusion coefficient should agree with Eq.( H2i) sufficiently before and after ts- 
Near ts the diffusion coefficient cannot be constant and the variance in position cannot be 
linear in time. Also the power is non-zero due to the induced electric field and the variance 
in velocity will be driven away from kT/m. 

In Figure [3] we show the numerically integrated thermodynamic quantities using the 
approximate Green's functions. Note the transient oscillations in the diffusion coefficient 
near ts and the expected agreement with Eq. (H2l) . For the protocol considered the change 
in field is oj {t) = asech^[a{t - ts)]. The power is positive when cr > and the work done 
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by the induced field causes the variance in velocity to spike above the equipartition value. 
Moreover, increasing the rate of the switch a induces a stronger electric field, resulting in 
a larger spike. Thus a slower rate corresponds to a smaller spike, which agrees with the 
expectation that in the adiabatic limit the system would be kept at equilibrium. When 
(T < the power is negative and the variance in velocity dips below the equipartition value. 

Since tanh is an odd function, the protocols for a forward = 3/2 and ahackward = — 3/2 are 
time-reversed with respect to each other. While the initial and final values of the diffusion 
coefficients match accordingly, it is evident from the figure that the transient process is not 
the same for the forward and backward protocols. There is a difference in the average power 
which is in fact the basis of testing non-equilibrium work theorems as in Ref . jl| . 



VII. CONCLUSION 



A time-dependent Langevin system may exhibit a stationary or non-stationary solution 
depending on the driving provided by the thermodynamic power P{t) = ^g{t)D{t). Equi- 
libration occurs when the power vanishes, implying that the variance in velocity takes the 
equipartition value. This happens if g{t) goes to zero at large times as in the case of the 
magnetic field switch, or if the diffusion coefficient goes to zero as in the case of a static linear 
potential. Non-equilibration happens when the power does not vanish in the long-time limit. 
The power could oscillate as in the case of a linear potential with a perturbative oscillating 
component, resulting in to thermodynamic quantities that oscillate as well. The power could 
also increase indefinitely, for example when g{t) = g < 0, resulting to unbounded variances. 
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FIG. 1: Thermodynamic quantities for an oscillating spring constant. The left (right) panel shows 
the variance in position (velocity) with the diffusion coefficient (power) in the inset. The black 
dashed line in the left (right) panel is the expected equilibrium value of (Ar^) = 1/g ((Aw?) = 
kT/m). 




FIG. 2: Comparison of the perturbative solution (dashed color) with the exact solution (solid 
color). The left (right) panel shows the variance in position (velocity) with the diffusion coefficient 
(power) in the inset. The black dashed lines in the diffusion coefficient inset are the values computed 
from Eq. (j42p . The black dashed line in the right panel is the free particle result, highlighting the 
exp[— 2/3t] decay. 
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FIG. 3: Thermodynamic quantities for a switched magnetic field. The left (right) panel shows the 
variance in position (velocity) with the diffusion coefficient (power) in the inset. The black dashed 
line in the diffusion coefficient inset is computed from Eq. ()42p . with the substitution to uj{t). 
The black dashed line in the right panel is the equipartition value kT/m. 



[1] A. Saha and A. M. Jayannavar, Phys. Rev. E. 77, 022105 (2008). 
[2] A. M. Jayannavar and M. Sahoo, Phys. Rev. E. 75, 032102 (2007). 
[3] M. Sahoo and A. M. Jayannavar, arXiv:0905.3901vl (2009). 
[4] F. Lillo and R. Mantegna, Phys. Rev. E. 61, R4675-R4678 (2000). 

[5] R. L. Liboff, Phys. Rev. 141, 222 (1966). This paper treats the equivalent system of a radial 

linear electric field crossed with an axial static magnetic field. 
[6] J. Dettman, Introduction to Linear Algebra and Differential Equations (Dover Publications, 

1986). 

[7] R. Mazo, in Stochastic Processes in Nonequilibrium Systems, edited by L. Garrido, P. Seglar 

and P. J. Shepherd (Springer- Ver lag, Berlin, 1978). 
[8] F. Paraan, M. Solon and J. P. Esguerra, Phys. Rev. E. 77, 022101 (2008). 



13 



